Crop Science

CoverCress is a new winter annual oilseed crop developed from field pennycress within the past 20 years. Field pennycress is commonly considered to be self-pollinated but little basic research has been published and there is some misalignment of conclusions. Our experience working with pennycress plant growth in greenhouse and field conditions over the past 13 years suggests that outcrossing is uncommon. We conducted lab, greenhouse, and field experiments to strengthen the body of work. Pollen viability kinetics analysis showed that longevity of pollen viability is negatively impacted by increasing temperatures and by direct exposure to light. Samples treated at 4C declined to 50% viability in 12 hours while it took just 2.5 hrs at 37C, and 1.6 hrs in full sunlight on a cool early April day. Cross-pollination was absent among greenhouse-grown plants flowering inside an agitated plastic pollen-containment covering. Across greenhouse tests, high rates of cross-pollination occurred only in an emasculation treatment that rendered flowers male sterile and opened the pistil to cross-fertilization. Field trials designed to measure pollen flow distance using a trackable fae1 knockout reporter gene failed to show detectable movement of pollen under field conditions in two locations. This data strongly suggests that domesticated field pennycress may be considered a self-pollinated crop and managed as such.

Spinach (Spinacia oleraceae) is a principal vegetable crop commercially grown in Controlled Environment Agriculture (CEA). Recent research suggests that root morphological and architectural differences among crop species influence yield, resource use efficiency, and environmental stress tolerance. These root traits may be exploited to increase yield, promote efficient nutrient use, and mitigate environmental stressors. This study measured differences between various spinach cultivars in CEA systems to reveal morphological and anatomical variation. We grew three spinach cultivars with different reported growing rates ( Income, Darkside, and El-Majestic) under NFT hydroponic and substrate-based systems in a controlled greenhouse environment over 45 days with destructive harvests at days 15, 30, and 45. Supplemental light (250 {micro}mol/m2/s) with 12-hour photoperiod and periodic fertigation was used. Harvests included the collection of leaf and root biomass, and scanning of root systems in WinRhizo software, measuring ten variables. On day 45, root cross-sections from orders 1-5 were embedded in JB-4 resin, sectioned, stained, and analyzed for diameter, vasculature, and rhizodermis characteristics. Results indicate that in spinach, differences in root system morphology are linked to cultivation systems over cultivar identity. Vascular and root anatomical alterations are minor compared to morphological differences in response to the cultivation system. Hydroponic-style growth systems are associated with the proliferation of fine-root ideotypes compared with substrate-based conditions. Such findings affirm previous studies, which suggest plastic root morphology in response to growth systems, and may be used to help create more resilient, resource-efficient cultivars. HighlightsO_LIIn spinach, root system morphology differences are linked to cultivation systems. C_LIO_LIRoot vascular and anatomical alterations are minor in response to cultivation system. C_LIO_LIHydroponic growth systems are linked to fine-root ideotype proliferation in spinach. C_LIO_LIFine-root ideotype proliferation may be a breeding target for CEA spinach. C_LI

Heat-tolerant rice cultivars are essential for mitigating global warming impacts. Basal anther dehiscence length (BDL) is a promising visible morphological marker for heat tolerance through stable pollination. We investigated the effects of sowing date on anther morphology, pollination, and fertility under controlled high-temperature conditions (35, 37, or 39 {degrees}C at flowering). Three japonica cultivars-- Akitakomachi (early heading), Koshihikari (medium), and Hatsushimo (late)--were sown monthly over 3 months and grown in pots. At heading, the plants were exposed to the temperature treatments for 3 days, and the proportion of florets with [≥]10 germinated pollen grains on the stigma (GP10) and seed set were assessed. Among anther traits, BDL showed the greatest variation, with all cultivars from the second sowing exhibiting the shortest BDL. Analysis of variance revealed significant effects of genotype, sowing date, and their interaction on anther traits and fertility. Regression analysis indicated that fertility was associated with GP10, with BDL contributing significantly to GP10 in the late-heading Hatsushimo, together with maximum temperature at flowering. Thus, both genotype and environment shape anther morphology, pollination, and fertility, indicating that BDL plasticity and genotype-specific environmental responses must be carefully considered when using BDL as a breeding marker for heat tolerance. HighlightVariation in sowing date significantly affects anther morphology and heat tolerance in rice. Genotype-specific responses to the growing environment require careful consideration for reliable breeding assessments.

Lodging in sorghum presents a significant challenge for plant breeders due to the trade-off between lodging resistance and grain yield. Manually measuring lodging across thousands of plots is time-consuming, expensive, and error-prone, making selection for lodging resistance challenging in breeding programs. Unmanned Aerial Vehicle (UAV) derived metrics offer a potential high-throughput, cost-effective alternative for lodging phenotyping. This study developed a framework for predicting plot-level lodging from UAV imagery across 2,675 sorghum breeding plots. Multi-temporal canopy height data were collected at two critical time points: maximum crop height and at manual lodging assessment. Height percentiles were extracted from UAV derived point clouds generated using photogrammetric algorithms. These data were used to develop parametric, non-parametric, and ensemble prediction models, which were evaluated using three statistical metrics. The ensemble model, averaging predictions from all models, achieved the highest accuracy with Pearson correlations of r = 0.80-0.84 and lowest residual mean square error (RMSE=16-18), explaining 64-70% of variation in manual lodging counts. Model diagnostics and iterative refinement, including inspection of UAV imagery and dataset curation, had minimal impact on model performance, demonstrating the robustness of the approach. Model performance was consistent across sites, with minimal effects of stratified sampling on accuracy, confirming the ensemble approach as optimal for plot-level lodging assessment. This study demonstrates that integrated multi-temporal UAV imagery offers a practical alternative to labor-intensive manual evaluation methods by enabling high-throughput lodging assessment suitable for implementation in sorghum breeding programs.

Improving photosynthesis is a promising approach to enhance sugar beet productivity. However, genetic variation in leaf photosynthesis and its relationship with disease resistance remain underexplored. We evaluated 98 sugar beet genotypes representing different breeding categories, including commercial F1 hybrids, seed-parent lines, and pollinator lines, in Hokkaido, northern Japan. Leaf gas exchange was measured during early growth under field conditions around the infection period of Cercospora leaf spot (CLS). To account for fluctuating irradiance during large-scale phenotyping, we applied a multilevel mixed-effects light-response model to estimate genotype-specific photosynthetic characteristics. Substantial genotypic variations in photosynthetic characteristics were detected. F1 hybrids exhibited higher photosynthetic capacity than breeding lines, whereas differences among breeding categories were unclear due to large within-category variation. Some breeding lines exhibited photosynthetic rates higher than those of hybrids, indicating exploitable genetic resources within the present genetic panel. We did not detect statistically significant trade-off between leaf photosynthesis and CLS resistance among 98 genotypes; in a subset of 19 genotypes analysed in detail, the relationship was even synergistic. Our results highlight the genetic diversity of leaf photosynthesis and its category-dependent structure, and suggest that selection for enhanced photosynthesis can proceed without substantial trade-off with CLS resistance. HighlightLeaf photosynthesis of 98 sugar beet genotypes showed significant genetic variation and dependence on breeding category. Active photosynthesis incurred minimal trade-off with Cercospora leaf spot resistance.

Phenotyping high-biomass perennial crops is laborious and the rate of genetic gain in perennial crop breeding programs is typically low. So, it is especially important to identify methods that produce efficiency gains in the breeding process. Miscanthus is a C4 perennial grass with favorable characteristics for producing biomass as a feedstock for biofuels and diverse biobased products. Increasing biomass yield will increase profitability and environmental benefits, so is a key target for Miscanthus breeding. In addition, the identification of well-adapted genotypes across a wide range of environmental conditions requires the establishment of multi-environment trials (METs). Sparse testing is a genomic prediction-based strategy that reduces the phenotyping costs in METs by selecting a subset of genotypes to evaluate in a subset of environments and then predicts the performance of the unobserved genotype-environment combinations. A Miscanthus sacchariflorus (MSA) population comprising 336 genotypes observed across three environments was analyzed. Three prediction models considering main effects (environments, genotypes, genomic) and interaction effects (genotype-by-environment; GxE interaction) were implemented for forecasting dry biomass yield (YDY), total culm (TCM), average internode length (AIL), and culm node number (CNN). Multiple calibration sets based on different compositions and sizes were considered to evaluate performance in terms of the predictive ability (PA) and the mean square error (MSE) for a fixed testing set size. The training set size ranged from 52 to 112 to predict a fixed set of 224 unobserved genotypes across all three environments. The results showed that the model accounting for GxE interaction presented the highest PA and the lowest MSE for CNN (PA: [~]0.77, MSE: [~]0.5) and YDY (PA: [~]0.70, MSE: [~]1.3) while for TCM and AIL these ranged from [~]0.28 to 0.41 and [~]1.3 to 4.3, respectively. Overall, varying training sets and allocation strategies did not affect PA and MSE, with 52 non-overlapping and 0 overlapping genotypes per environment as the optimal cost-effective allocation framework. This suggests that implementing sparse testing designs could significantly reduce phenotyping costs by fivefold, without compromising PA in breeding programs for perennial crops such as Miscanthus.

The genetics of interspecific groups remains largely unexplored, despite the central role of social (or indirect) genetic effects in shaping phenotypic expression within communities. Intercropping, i.e. the simultaneous cultivation of multiple crop species in the same field, offers a powerful model to harness these interspecific social effects. Such species mixtures provide well-documented agricultural benefits, yet few breeding frameworks have integrated the genetics of social interactions. Here, we address this gap by extending quantitative genetic theory to interspecific groups, with intercropping as a concrete and applied model case. We propose a quantitative genetic model that jointly analyzes intra and interspecific interactions within a unifying framework. Breeding values are decomposed into a direct component, shared in mono and mixed-crops, an interspecific social component corresponding to the effect of one species on another, and an intraspecific component that captures the social effects within a mono-genotypic stand of cloned plants. Statistically, this consists in simultaneously fitting several linear mixed models, one per stand type, all having direct breeding values in common. As no open-source software can fit such a complex mixed model, we provide such an implementation in R/C++. Simulations across various genetic (co)variance structures and sparse experimental designs showed accurate estimation of all genetic (co)variances and breeding values. With an incomplete, yet balanced design combining sole crops and intercrops, genetic gains in both systems were achievable simultaneously, enabling breeding strategies that progressively integrate intercropping into existing, sole-crop-only schemes. More broadly, this framework allows dissecting direct and social genetic effects when genotypes are observed in mono- and mixed-species situations, cultivated or not.

Drought is a significant factor in agricultural losses, making it imperative to understand how root system architecture (RSA) adapts to environmental condition like water deficit. HydroRoot is a functional-structural plant model (FSPM) aimed at analyzing and simulating hydraulic and solute transport of RSA. The model integrates a static hydraulic solver, a coupled water-solute transport solver, a statistical generator of RSA based on Markov model, and a dynamic hydraulic model accounting for root growth. This paper presents the model, the mathematical description of the formalism of solvers, and use cases with their associated tutorials. Five use cases illustrate capabilities of HydroRoot, which has been successfully used for phenotyping root hydraulics across various species, including Arabidopsis, maize, and millet. The model-driven phenotyping method "cut and flow" is presented to characterize axial and radial conductivities on a given root genotype. Finally, three step-by-step tutorials provide a structured way to learn how to use HydroRoot 1) to simulate hydraulic on a given architecture, 2) to simulate water and solute transport on a maize root, and 3) to simulate hydraulic on two pearl millet genotypes with varying soil conditions. Hydroroot is an open-source package of the OpenAlea platform, with the code publicly available on Github. A comprehensive documentation is available with a reproducible gallery of examples.

Global networks of crop breeding programs leverage diverse germplasm, but diversity increases the complexity of maintaining stability in their elite genepools. To characterize genetic heterogeneity in breeding metapopulations and develop insights on how to manage it, we simulated the evolution of breeding populations on fitness landscapes. We revealed the geometric decrease in the average effect size of alleles segregating as standing variation that become fixed along an adaptive walk. We also demonstrated how independent adaptive walks of subpopulations are influenced by genetic drift, leading to cryptic genetic heterogeneity among elite genepools. This variation is released when elite lines derived from independent subpopulations are crossed, leading to segregation for 2-4X more major QTL in admixed families as in unadmixed families, and 2-4X more epistatic interactions. The emergent property of fitness epistasis for traits under stabilizing selection is well-understood in evolutionary genetics, but under-appreciated in crop quantitative genetics. To highlight the importance of this phenomenon, we constructed an empirical genotype-to-fitness landscape from the sorghum NAM, a global admixed prebreeding resource, demonstrating the utility of fitness landscapes for inferring genetic compatibilities within metapopulations. Our findings suggest that in breeding networks, strategies for effective germplasm exchange must account for epistasis in the oligogenic component of the genetic architecture of locally-adapted traits. Article summaryModern public sector crop improvement happens in networks of breeding programs that routinely exchange genetic information. Traditional models for understanding quantitative traits have limited predictiveness in situations with such genetic heterogeneity. This study uses breeding simulations and empirical data to show the utility of the fitness landscape framework for characterizing the genetic architecture of complex traits in breeding metapopulations. By simulating the evolution of breeding programs and integration into networks, it demonstrates how epistatic interactions between large-effect alleles are a fundamental property that must be accounted for when exchanging germplasm. Graphical Abstract O_FIG O_LINKSMALLFIG WIDTH=200 HEIGHT=102 SRC="FIGDIR/small/712732v1_ufig1.gif" ALT="Figure 1"> View larger version (25K): org.highwire.dtl.DTLVardef@1541326org.highwire.dtl.DTLVardef@b553a8org.highwire.dtl.DTLVardef@8758b4org.highwire.dtl.DTLVardef@1d0bdcd_HPS_FORMAT_FIGEXP M_FIG C_FIG

Grapevine cluster architecture is a key selection target in breeding programs because it influences disease susceptibility, yield stability and juice quality. High-throughput phenotyping offers a rapid and non-destructive approach to capture biochemical and structural variation in these traits, yet the influence of plant organ reflectance and data partitioning strategies on trait prediction remains poorly understood. In this study, we evaluated how hyperspectral reflectance from different grapevine organs contributes to the prediction of cluster architecture and juice quality traits in two clonal populations of Riesling and Pinot. Using partial least squares regression (PLSR), we assessed the prediction accuracy of eight cluster architecture and six juice quality traits under two data partitioning strategies. Models based on cluster reflectance outperformed those using dry leaf reflectance for most traits, except for pH. Partitioning the dataset by cluster type increased trait variance and improved predictions for number of berries (R{superscript 2} = 0.53), berry diameter (R{superscript 2} = 0.79), and total acidity (R{superscript 2} = 0.48). Visible, red-edge and NIR spectra were most informative regions to predict the traits studied. Together, our results highlight the importance of organ-specific data and appropriate calibration strategies to improve phenomic models for the development of scalable proxies for grapevine improvement. HighlightSpectral phenomics reveals that prediction accuracy in grapevine depends on organ spectral signatures and traits, with cluster reflectance outperforming leaves, informing new phenotyping strategies for breeding improvement.

Genomic selection holds the potential to serve as a strategic tool to enhance the genetic gain of complex traits in Miscanthus breeding programs. The development of improved cultivars requires their assessment for various traits across diverse environments to ensure suitable overall performance. Hence, the multi-trait multi-environment (MTME) genomic prediction (GP) models offer an opportunity to improve selection accuracy. This study aims to evaluate the potential of five GP models: (1) three MTME models including genotype-by-trait-by-environment interaction (GxExT) and (2) two single-trait multi-environment (STME) models (with and without GxE interaction). A Miscanthus sacchariflorus population comprising 336 genotypes evaluated in three environments and scored for four traits (biomass yield YDY, total culm number TCM, average internode length AIL, and culm node number CNN) was analyzed. The predictive ability of the models was evaluated considering three cross-validation schemes resembling realistic scenarios (CV1: predicting new genotypes, CVP: predicting missing traits in a given environment, and CV2: predicting partially observed genotypes). On average, in all cross-validation schemes compared to the STME the predictive ability of the MTME models was 10% to 70% higher for TCM and AIL. On the other hand, for YDY and CNN, both STME models performed similarly or slightly better (between 5 to 64%) than the MTME models in most environments. While the MTME models were not successful for all traits when compared to their STME counterparts, MTME models improved the prediction of the performance of genotypes that were untested across environments or lacked trait information in a specific environment. Overall, our study suggests that MTME GP models can be implemented in Miscanthus breeding programs to improve the predictive ability of the complex traits, shorten breeding cycles, and accelerate selection decisions.

Potato common scab (Streptomyces sp.) is an economically important disease that reduces the quality and market value of tubers. A key aspect in developing management strategies involves accurately quantifying the disease. Due to the three-dimensional nature of the tuber and the heterogeneous distribution of lesions across its surface, visual estimates of severity can be challenging. Therefore, the objectives of this study were to develop and validate a standard area diagram (SAD) for estimating common scab severity on potato tubers and to compare validation outcomes obtained using real tubers and digital images. A SAD comprising six severity levels (from 1.3 to 66.8%) was developed based on image analysis of naturally infected tubers. Validation was conducted using two complementary approaches in which inexperienced raters evaluated either real potato tubers or digital images of the same tubers under unaided and aided conditions. Accuracy, bias components, and inter-rater reliability were quantified using absolute error metrics, Lins concordance correlation coefficient, intraclass correlation coefficients, and overall concordance correlation coefficients. Use of the SAD significantly improved accuracy, reduced systematic bias, and increased inter-rater reliability across both validation approaches. No significant differences were detected between assessments conducted on real tubers and images, although image-based evaluations showed a slight, non-significant tendency toward reduced scale and location bias under aided conditions. These results demonstrate that a dimension-aware SAD integrating information across the full tuber surface enhances the reliability and reproducibility of visual severity assessments and supports the use of image-based evaluations for training, large-scale surveys, and remote or collaborative applications involving three-dimensional plant organs.

This study investigated the mutagenic effects of ethyl methane sulfonate (EMS) on the M{square} generation in cowpea (Vigna unguiculata (L.) Walp.) cultivar Wang Kae. A total of 275 M{square} seeds were treated with EMS concentrations of 20 mM, 40 mM, and 80 mM (75 seeds per treatment) by soaking for six hours, while 50 untreated seeds served as the control (0 mM). Phenological, yield-related and yield traits were recorded, and data were analysed using Jamovi 2.7.15 and JASP 0.95.4.0 through one-way ANOVA with post hoc contrast, principal component biplot, and cluster analyses. No optimal mutagenic concentration (LD50) was identified. Seed germination and seedling survival rates increased with increasing EMS concentration, ranging from 70.00% and 62.00% in the control (0 mM) to 89.33% and 74.67% at 80 mM, following the trend 0 mM < 20 mM < 40 mM < 80 mM. Significant differences (P < 0.05) were observed among treatments for all phenological traits, pod length, locule number, seed traits, and yield per plant. Yield was significantly higher (P = 0.047) at 20 mM (61.19 {+/-} 3.34 g) compared to the control. Contrast analysis identified genotypes B33 and D56 as the most productive mutants, with yields of 125.44 g and 111.85 g, respectively. Principal component analysis extracted eighteen components, with the first four cumulatively explaining 50.60% of total variation. Biplot analysis of PC1 and PC2 captured all phenological traits, key seed traits, and yield attributes, highlighting the superior performance of B33 and D56. Cluster analysis partitioned the 190 genotypes into six groups, with B33 and D56 constituting distinct clusters. EMS mutagenesis effectively induced heritable phenotypic variation, with putative superior genotypes identified for advancement to M{square} and evaluation in replicated multi-environment trials toward the development of farmer- and consumer-preferred cowpea varieties.

In cereal-legume intercrops, weed suppression is primarily driven by cereals, whose competitiveness is shaped by trait plasticity--morphological adjustments in response to the intercrop environment. However, how individual cereal traits respond plastically and contribute to system performance remains unclear, hampering improvements through breeding or system design. We combined field experiments with functional-structural plant modelling to quantify plastic responses of four cereal traits (tiller number, tiller angle, specific leaf area (SLA), and specific internode length (SIL)) and their effects on weed suppression and crop productivity. Field measurements revealed plasticity in tiller number, tiller angle, and SIL between sole crops and intercrops, while SLA showed minimal differences. Simulations showed that intermediate tiller numbers resulted in the strongest weed suppression and highest productivity, indicating an optimum, while more horizontal tillers suppressed weeds slightly better than vertical ones. Weed suppression increased with higher SLA values, while SIL showed a saturating response, increasing to intermediate SIL values and plateauing thereafter. In simulations with short-statured cereal phenotypes (low SIL), the reduction in cereal weed suppression was compensated by the legume component. This study demonstrates how FSP modelling can be used to investigate trait plasticity mechanisms and generate testable hypotheses about trait effects in complex intercrop systems. HighlightCereal trait plasticity shapes weed suppression in cereal-legume intercrops, with distinct response patterns per trait, while legumes can compensate for weakly competitive cereals, suggesting balanced competition over cereal dominance.

Bioavailability of iron, an essential micronutrient to plants, is low in alkaline or calcareous soils, which are prevalent across semi-arid production regions. Breeding efforts to increase tolerance to iron deficiency chlorosis (IDC) in sorghum, a major crop of semi-arid regions, are confounded by spatial variation of stress severity in field trials. Here we developed and validated two high-throughput phenotyping approaches to address this challenge, with multi-spectral aerial imaging in the field and a controlled-environment assay to isolate the effects of iron bioavailability. In the field, severity and uniformity of stress are highly predictive of genetic signals for IDC tolerance (R2 > 0.6 for soil pH metrics and H2). Plot-level data filtering for stress conditions based on control genotypes successfully addresses field spatial variation (unfiltered H2 = 0.18 vs. filtered H2 = 0.4). The controlled-environment assay proxies field stress using iron sources with differential bioavailability, evidenced by high heritability ( H2 = 0.98) and phenotypic differential for hybrid control genotypes that matches field performance. Finally, we show that assay phenotypes are suitable for genome-wide association studies in global germplasm. Together, these field and lab phenomic approaches can be deployed to understand genetics of IDC tolerance and develop crops resilient to alkaline soils. HIGHLIGHTStress severity and uniformity greatly impact detection of genetic signals underlying iron deficiency chlorosis tolerance in sorghum. A controlled-environment assay reduces spatial heterogeneity and improves assessment of tolerance genetics.

Enset (Ensete ventricosum), a multipurpose crop domesticated exclusively in Ethiopia, serves as a staple food for millions of smallholder farmers. It is primarily cultivated as a monocrop in homegardens, leaving it vulnerable to climate change risks. One potential nature-based solution involves agroforestry systems; however, ensets response to canopy cover remains unclear. This study examined how scattered trees in enset farms affected microclimate and enset morpho-physiology in South Ethiopia. Trees significantly modified the microclimate conditions in enset homegardens. The average daily reductions in air, soil surface, and soil temperatures ranged from -0.5 to -1.9 {degrees}C, -0.4 to -2.1 {degrees}C, and +0.4 to -1.0 {degrees}C, respectively. The minimum soil moisture offset ranged from +0.8% to +5.7%. Although the tree identity effect on enset growth was negligible, planting position relative to the overstory trees significantly influenced enset responses. Most morphophysiological traits were higher under tree canopies, with progressively lower values at the edge and outside the tree canopy. In contrast, leaf dry matter content exhibited an inverse trend, aligning with the leaf economics spectrum. These results demonstrate ensets adaptability to canopy shade, suggesting potential for agroforestry expansion. Cultivar-specific shade tolerance and ideal shade levels to maintain enset productivity should be investigated further.

Deciphering the genetic architecture of complex quantitative phenotypes remains challenging in quantitative genetics. These traits not only depend of multiple genetic factors but are also established over time and environments. Although quantitative genetics has investigated the genetic determinism of phenotypic plasticity in contrasted environmental conditions, the time related phenotypic plasticity has received less attention. Here we proposed a multivariate Bayesian framework, the Bayesian Varying Coefficient Model, designed for analysing the genetic architecture of the time related phenotypic plasticity by a multilocus approach. We applied the BVCM to time series phenotypes measured at various time scales (daily, monthly, yearly) across a diverse set of biological species. We included in this study: yeast (Saccharomyces cerevisiae), fungi (Fusarium graminearum), eucalyptus (Eucalyptus urophylla x E. grandis), and sweet cherry tree (Prunus avium). The BVCM results were compared with those obtained with a known genome-wide association method carried out time by time. For all species and traits, the BVCM was able to detect the major QTL identified by marker-trait association methods and revealed additional genetic regions of weak effect. It also increased the phenotypic variance explained for most of the phenotypes considered. It revealed dynamic QTLs with transitory, increasing or decreasing effects over time. By considering both the temporal and genetic multivariate structures in a single statistical model, we increased our understanding of the genetic architecture of complex traits notably by reducing the issue of missing heritability. More broadly, this work raises the foundation for extended applications in functional genomics, evolutionary ecology, and crop breeding programs, in which time-related phenotypic plasticity remains crucial for predicting and selecting key quantitative complex traits. Key messageBy capturing the genetic factors influencing the time related phenotypic plasticity, our approach contributes to a deeper understanding of the dynamic nature of genotype-phenotype relationships.

Understanding crop genetic diversity is essential for conservation and breeding, yet farmer-maintained germplasm remains largely underrepresented in genomic studies. Theobroma cacao L. has a complex domestication history and extensive global diversity, and cacao currently cultivated in Central America, particularly in Costa Rica, has been understudied compared to South American and Mexican cultivars despite cultural and historical importance. In this study, we investigate the genetic diversity of cacao from farmer-managed systems across Costa Rica to search for Criollo germplasm and identify and characterize any unique local genetic groups. Ninety-four trees were sampled from 17 farms across four regions of the country and sequenced using whole genome resequencing. Farmer materials were analyzed alongside 166 previously characterized reference accessions representing major cacao genetic groups. Population structure analyses, phylogenetic reconstruction, and network approaches revealed that Costa Rican cacao encompasses multiple known genetic groups, including Criollo-derived lineages, while also harboring locally distinct diversity not fully represented in current global reference collections. Analyses revealed close kinship between many accessions with no clear geographic patterns corresponding to the observed population differentiation, reflecting the effects of farmers in creating dominant patterns of gene flow through seed-saving, clonal propagation, and sharing genotypes among farms. Heterozygosity levels varied substantially among individuals, consistent with a mixture of highly inbred Criollo trees and more heterozygous, admixed genotypes. We find that farmer-managed cacao systems are reservoirs of genetic diversity, including possibly rare or historically important lineages, underscoring the value of these farming systems for effective conservation and management of genomic resources for cacao resilience and improvement.

Salinity and sodicity stresses adversely affect rice growth and yield. To overcome yield losses, suitable tolerant rice cultivars can be developed through a marker-assisted breeding (MAB) program. In the present study, genomic regions associated with sodicity stress tolerance at the reproductive stage were identified using a high-density 50kSNP array in a recombinant inbred line (RIL) population derived from the contrasting rice genotypes CSR11 and MI48. A total of 50 QTLs were detected for various yield-related traits; further, 19 QTLs with [&ge;]15% of phenotypic variance were selected for integrated (omics) analysis. RNA sequencing of leaves and panicles at the reproductive stage under sodic stress conditions was employed to find differentially expressed genes. A total of 1368 and 1410 SNPs; 104 and 144 indels were found for MI48 and CSR11, respectively, within the QTL regions from resequencing. At chromosomes 1 and 6, colocalized QTLs (qPH1-1/qGP1-1 and qGP6-2/qSSI6-2) were discovered. Differentially expressed genes (DEGs) were mapped over the QTL regions selected, and SNP variations and indels were screened for colocalized QTLs. Potential candidate genes, namely Os-pGlcT1 (Os01g0133400), OsHKT2;1 (Os06g0701600) and OsHKT2;4 (Os06g0701700), OsANTH12 (Os06g0699800), and OsPTR2 (Os06g0706400), were identified as being responsible for glucose transport, ion homeostasis, pollen germination, and nitrogen use efficiency, respectively, under salt stress. Finally, our study provides important insights into the genes and potential mechanisms affecting grain yield under sodic stress in rice, which will contribute to the development of molecular markers for rice breeding programs.

Black spruce (Picea mariana [Mill.] B.S.P.) is an emblematic and ubiquitous species of the North Americas boreal forest. While conifer breeding programs have traditionally focused on growth and wood property traits, the study of climate adaptation traits is becoming increasingly prevalent, given the predicted impact of climate change on North Americas boreal zone. Through this study, we aimed to identify genes associated with climate adaptation in black spruce across Canada. A total of 254 black spruce trees from 30 populations, covering most of the species distribution range, were sampled and genotyped for SNPs located in [~]5000 gene loci. Uni- and multivariate Genotype-Environment Association (GEA) approaches, namely LFMM and RDA, as well as an outlier method based on population differentiation (FST) were used to identify genes significantly associated with climatic factors. As such, a total of 77 genes carrying significant candidate SNPs were identified, among which 14 candidates were corroborated by at least two methods. Many of these gene SNPs were also confirmed at a smaller geographic scale, across west - east partitions corresponding to the two main black spruce historical lineages. Notably, significant gene SNPs were more frequently associated to moisture/aridity factors in the western part of the range, and more to temperature factors in the eastern part. The genes carrying these SNPs were also frequently associated to abiotic and biotic stress response. In the context of rapid climate change in the Canadian boreal forest, the results obtained within the framework of this study should support implementing gene conservation efforts while assisting prediction in black spruce breeding programs, which are instrumental to producing adapted planting stock for the large-scale reforestation efforts conducted annually across the Canadian boreal forest.